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In this work we analyse the possibility of constraining the QCD dynamics at high energies studying the heavy 
quark photoproduction at LHC in coherent interactions. The rapidity distribution and total cross section for 
charm and bottom production are estimated using three different phenomenological saturation models which 
successfully describe the HERA data. Our results indicate that the experimental study of the inclusive heavy 
quark photoproduction can be very useful to discriminate between the classical and quantum versions of the Color 
Glass Condensate (CGC) formalism. 
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1. INTRODUCTION 

The cross sections for heavy quark production 
in hadron-hadron and lepton-hadron collisions at 
high energies are strongly dependent on the be- 
havior of the gluon distribution, which is deter- 
mined by the underlying QCD dynamics. The- 
oretically, at high energies the QCD evolution 
leads to a system with high gluon density, charac- 
terized by the limitation on the maximum phase- 
space parton density that can be reached in the 
hadron wave function (parton saturation). The 
transition is specified by a typical scale, which 
is energy dependent and is called the satura- 
tion scale. Signals of parton saturation have al- 
ready been observed both in ep deep inelastic 
scattering at HERA and in deuteron-gold colli- 
sions at RHIC. However, the observation of this 
new regime still needs confirmation and so there 
is an active search for new experimental signa- 
tures. In this contribution we study the in- 
clusive and diffractive photoproduction of heavy 
quarks in coherent proton-proton collisions con- 
sidering three phenomenological models based on 
the Color Glass Condensate formalism, which de- 
scribe quite well the current experimental HERA 
data for inclusive and exclusive observables (For 
a detailed discussion see Ref. pQ). 

In hadron-hadron colliders, the relativistic pro- 



tons give rise to strong electromagnetic fields, 
which can interact with each other. Namely, 
quasi-real photons scatters off protons at very 
high energies in the current hadron colliders [3J. 
Our main motivation comes from the fact that in 
coherent interactions at Tevatron and LHC the 
photon reaches energies higher than those cur- 
rently accessible at DESY - HERA. The heavy 
quark photoproduction cross section in a proton- 
proton collision is given by, 



a(p + p^p + QQ + Y) = 2 
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where u is the photon energy in the center-of- 
mass frame (c.m.s.), W- 



ii' 



is the c.m.s. photon- 
proton energy and \J Snn denotes the proton- 
proton c.m.s. energy. The final state Y can be a 
hadronic state generated by the fragmentation of 
one of the colliding protons (inclusive production) 
or a proton (diffractive production). 

2. PHOTOPRODUCTION OF HEAVY 
QUARKS 

In the color dipole approach the inclusive and 
diffractive heavy quark photoproduction cross 
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section are given by [11314) 

atot =2 Jd 2 b (Af(x,r,b)) (2) 
and 

ag t = J d 2 b {M 2 (x,r,b)) . (3) 
where 

<(...)) = J d 2 r J dz|vl/ 7 (r,z)| 2 (...) (4) 

In these equations the variable r defines the rel- 
ative transverse separation of the pair (dipole), 
z (1 — z) is the longitudinal momentum fractions 
of the quark (antiquark) and the function \l/ 7 is 
the light-cone wave function for transversely po- 
larized photons, which depends in our case of 
the charge and mass of the heavy quark. The 
function N(x,r,b) is the forward dipole-target 
scattering amplitude for a dipole with size r and 
impact parameter b which encodes all the infor- 
mation about the hadronic scattering, and thus 
about the non-linear and quantum effects in the 
hadron wave function [5]. During the last years 
an intense activity in the area resulted in sophisti- 
cated models for the dipole-proton scattering am- 
plitude, which have strong theoretical constraints 
and which are able to describe the HERA and/or 
RHIC data. In what follows we will use three dis- 
tinct phcnomenological saturation models based 
on the Color Glass Condensate formalism which 
describe quite well the more recent HERA data: 
the IIM [6], the bCGC [7J and the IP-SAT 
models. In the bCGC model, which is the quan- 
tum version of the CGC formalism, the dipole- 
proton scattering amplitude is given by 




Furthermore, we will use in our calculations the 
scattering amplitude scattering proposed in Ref. 
[5], denoted IP-SAT, which is given by 



JV^ AL (x,r,b) = 1-S(x,r,b) , (6) 



where 

S = cxp (--^-r 2 a s (fi 2 ) xg(x, fi 2 ) T(6)) , (7) 

the scale is related to the dipole size r by 
fi 2 = 4/r 2 + /Uq and the gluon density is evolved 
from a scale Mo up to /i 2 using LO DGLAP evo- 
lution without quarks assuming that the initial 
gluon density is given by xg(x, = A g x~ Xg (1— 
a;) 5 ' 6 . The values of the parameters /i 2 ,, A g 
and X g are determined from a fit to Fi data. 
Moreover, it is assumed that the proton shape 
function T(b) has a Gaussian form, T(b) = 
l/(27rB G )exp[-(b 2 /2B G )], with B G being a free 
parameter which is fixed by the fit to the differ- 
ential cross sections for exclusive vector meson 
production. The parameter set used in our cal- 
culations is the one presented in the first line of 
Table III of [Sj: Mo = 1-17 GeV 2 , A g = 2.55, 
A 9 = 0.020 and Bq — 4 GeV" 2 . The previous 
expression for the forward scattering amplitude 
can be obtained to leading logarithmic accuracy 
in the classical effective theory of the Color Class 
Condensate formalism. Moreover, it is applica- 
ble when the leading logarithms in Q 2 dominate 
the leading logarithms in 1/x, with the small r 
limit being described by the linear DGLAP evolu- 
tion at small- ir. In contrast, the bCGC model for 
M captures the basic properties of the quantum 
evolution in the CGC formalism, describing both 
the bremsstrahlung limit of linear small- a; evo- 
lution (BFKL equation) as well nonlinear renor- 
malization group at high parton densities (very 
small- re). Consequently, the IP-SAT model can 
be considered a phcnomenological model for the 
classical version of the CGC, while the bCGC for 
the quantum limit. It is important to emphasize 
that both models provide excellent fits to a wide 
range of HERA data for x < 0.01. Therefore, 
the study of observables which are strongly de- 
pendent on Af is very important to constrain the 
underlying QCD dynamics at high energies. In 
what follows we consider these two models as in- 
put in our calculations of the inclusive and diffrac- 
tive heavy quark photoproduction in coherent pp 
interactions at LHC energies. 
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Figure 1. The rapidity distribution for the in- 
clusive charm (upper panel) and bottom (lower 
panel) photoproduction on pp reactions at LHC 
energy ^/Snn = 14TeV. Different curves corre- 
spond to distinct phenomenological models. 



3. RESULTS 

The distribution on rapidity y of the produced 
open heavy quark state can be directly computed 
using its relation with the photon energy to, i.e, 
y oc ln-£-. A reflection around y = takes into 
account the interchange between the proton emit- 
ting the photon and the target proton. Explicitly, 
the rapidity distribution is written down as, 

da [pp^pQQY)] _ dN 7 (u) 

Ty ^ a iP~*QQ Y1 W 

where Y is a hadronic final state X produced by 
proton fragmentation in the inclusive case and 



Y = p for diffractive production. 

The resulting rapidity distributions for inclu- 
sive and diffractive heavy quark photoproduction 
at LHC energies coming out of the distinct phe- 
nomenological models considered in previous sec- 
tion are depicted in Figs. [T] and [2l respectively. 
For the inclusive case (Fig. [TJ, the IIM and 
bCGC predictions are very similar. In contrast, 
these predictions are distinct in the diffractive 
case (Fig. with the bCGC prediction being 
larger than IIM one at mid-rapidity. On the other 
hand, the IP-SAT prediction is larger than these 
predictions by a factor 2 (3) in the charm (bot- 
tom) case. We can consider the IIM and bCGC 
predictions as lower bounds for the coherent pro- 
duction of heavy quarks at LHC. Our results indi- 
cate that the experimental study of the inclusive 
heavy quark photoproduction can be very useful 
to discriminate between the classical and quan- 
tum versions of the CGC formalism. It also is 
true in the diffractive case, where the different 
models can be discriminated more easily. 

The results for the integrated cross section con- 
sidering the distinct phenomenological models are 
presented in Table [TJ for the inclusive and diffrac- 
tive charm and bottom pair production at LHC. 
The IP-SAT model gives the largest rates among 
the approaches studied, followed by the bCGC 
and IIM models with almost identical predictions. 
Therefore, these reactions can have high rates at 
the LHC kinematical regime. The cross sections 
for diffractive production are approximately two 
orders of magnitude smaller than the inclusive 
case, but due the clear experimental signature of 
this process (two rapidity gaps), its experimental 
analysis still is feasible. 

4. CONCLUSIONS 

In summary, we have computed the cross sec- 
tions for inclusive and diffractive photoproduc- 
tion of heavy quarks in pp collisions at LHC ener- 
gies. This has been performed using modern phe- 
nomenological models based on the Color Glass 
Condensate formalism, which describe quite well 
the inclusive and exclusive observables measured 
in ep collisions at HERA. The obtained values arc 
shown to be sizeable at Tevatron and are increas- 
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Table 1 

The integrated cross section for the inclusive and 
diffractive photoproduction of heavy quarks in pp 
collisions at LHC energies. 




Figure 2. The rapidity distribution for the 
diffractive charm (upper panel) and bottom 
(lower panel) photoproduction on pp reactions at 
LHC energy ^/Snn = 14TeV. Different curves 
correspond to distinct phenomenological models. 



ingly larger at LHC. The feasibility of detection 
of these reactions is encouraging, since their ex- 
perimental signature should be sufficiently clear. 
Furthermore, they enable to constrain the under- 
lying QCD dynamics at high energies, which is 
fundamental to predict the observables which will 
be measured in central hadron-hadron collisions 
at LHC [9]. 

REFERENCES 

1. V. P. Goncalves, M.V.T. Machado and A.R. 
Meneses, Phys. Rev. D. 80 (2009) 034021. 

2. G. Baur, K. Henckcn, D. Trautmann, S. 



Sadovsky, Y. Kharlov, Phys. Rep. 364 (2002) 
359; C. A. Bertulani, S. R. Klein and J. Nys- 
trand, Ann. Rev. Nucl. Part. Sci. 55 (2005) 
271. 

3. V. P. Goncalves and M.V.T. Machado, Phys. 
Rev. D 71 (2005) 014025. 

4. V. P. Goncalves and M.V.T. Machado, Phys. 
Rev. D 75 (2007) 031502. 

5. E. Iancu and R. Venugopalan, 
|arXiv:hep-ph/0303204| A. M. Stasto, 
Acta Phys. Polon, B 35 (2004) 3069; H. 
Weigert, Prog. Part. Nucl. Phys. 55 (2005) 
461; J. Jalilian-Marian and Y. V. Kovchegov, 
Prog. Part. Nucl. Phys. 56 (2006) 104. 

6. E. Iancu, K. Itakura, S. Munier. Phys. Lett. 
B 590 (2004) 199. 

7. G Watt and H. Kowalski, Phys. Rev. D 78 
(2008) 014016. 

8. H. Kowalski, L. Motyka and G. Watt, Phys. 
Rev. D 74 (2006) 074016. 

9. N. Armesto et al, J. Phys. G 35 (2008) 
054001. 



